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We study the quantum Rabi model within the framework of the analytical solution developed in Phys. Rev. 
Lett. 107, 100401 (201 1). In particular, through time-dependent correlation functions, we give a quantitative cri- 
terion for classifying two regions of the quantum Rabi model, involving the Jaynes-Cummings, the ultrastrong, 
and deep strong coupling regimes. In addition, we find a stationary qubit-field entangled basis that governs the 
whole dynamics as the coupling strength overcomes the mode frequency. 
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I. INTRODUCTION 

The light-matter interaction has been subject of central in- 
terest since the early years of quantum mechanics. In 1936, 
one of the first attempts to explain the results coming from 
experiments was the Rabi model, that describes the sim- 
plest dipole semiclassical interaction between light and mat- 
ter Q, and is reduced to a pseudospin-1/2 system driven by a 
monochromatic classical radiation field. However, the advent 
of quantum technologies such as cavity QED ElO, have al- 
lowed us to access the quantum regime of the radiation field, 
where the dynamical description is given by the celebrated 
Jaynes-Cummings model (JCM) [4]. This model predicts col- 
lapses and revivals of the population inversion, the appearance 
of Jaynes-Cummings doublets as a consequence of the excita- 
tion number conservation, and it has found a testbed in several 
hybrid setups such as trapped ions [5], quantum dots (6], and 
circuit quantum electrodynamics (circuit QED) I7UT611. 

Circuit QED has been growing both theoretically and ex- 
perimentally, and their complex proposals such as the gener- 
ation of multipartite entanglement fT7lfT9l rely on the fun- 
damentals of the JCM and the Tavis-Cummings model l20l . 
This can be justified because ratios between the coupling 
strength g and the resonator frequency u may grow from typ- 
ical quantum optical values of g/uo ~ 10 -6 to circuit QED 
values of g/oj ~ 10 -2 . Note that, in the latter, the rotating- 
wave approximation (RWA) can still be applied. Nowadays, 
two key experiments with superconducting circuits I2T1 l22l 
have made a significant improvement in the coupling strength, 
reaching values g/uj ~ 0.1 in the so-called ultrastrong cou- 
pling (USC) regime l23lJ25l In this case, the RWA is not 
longer valid and all dynamical and statical properties have to 
be explained through the quantum Rabi model (QRM) 

Hr = hAa z + ftwa^a + hga x (a + a^), (1) 

where a z and a x are Pauli matrices, a{a)) is the annihilation 
(creation) operator, A is the qubit energy, uo is the resonator 
frequency, and g stands for the coupling strength. In addi- 
tion, a recent proposal considers the case where the coupling 
strength g becomes comparable or larger than the mode fre- 
quency uj, g/u > 1, which is called deep strong coupling 
(DSC) regime [26]. In this case, the dynamics can be intu- 
itively explained as photon number wavepackets propagating 



along a defined parity chain. Although the state-of-the-art in 
circuit QED does not provide this coupling strength yet, the 
main features of the DSC regime have been observed in an 
analog quantum simulation I27ll28l . 

The QRM described by the Hamiltonian ([I]), has substan- 
tial differences as compared to the JCM, in fact only re- 
cently the properties of the QRM have been completely un- 
derstood (291I30). In the QRM, a discrete Z 2 -symmetry re- 
places the continuous U (l)-symmetry of the JCM. Therefore 
the excitation number is no longer a conserved quantity and 
the Hilbert space splits into two infinite-dimensional invariant 
subspaces, the parity chains [26]. Each eigenstate can be la- 
beled with a Z2 -quantum number, the parity P — a z e l7Taia , 
taking values ±1. 

Furthermore, it is possible to give analytical expressions for 
these eigenstates as elements of the Bargmann space [ 31 ], that 
allows the well-defined computation of their norms and over- 
laps with eigenstates of the harmonic oscillator, without trun- 
cation of the Hilbert space 1321 . 

The aim of this work is to present a new insight of the quan- 
tum Rabi model by stressing the use of dynamical correla- 
tion functions coming from the analytical solution obtained in 
Ref. [29]. In this manner, we are able to explain relevant fea- 
tures such as the validity of the RWA in two well-defined re- 
gions, ranging from the JC regime to higher-coupling regimes 
of the quantum Rabi model. In addition, as the coupling 
strength g enters into the DSC regime, we find that the true 
eigenstates of the system can be well approximated by the 
shifted oscillator basis in each invariant parity chain with- 
out the need for a more complicated basis. This fact is sup- 
ported by the calculation of the Wigner function of the eigen- 
states, whose unexpected fidelity can be understood in terms 
of the analytical form of the eigenfunctions |29lL resembling 
closely to Fock states. In this context, we also find station- 
ary Schrodinger cat-like states in the DSC regime. Finally, we 
present our concluding remarks. 



II. DYNAMICAL CLASSIFICATION 

In this section, we shall use the exact dynamics of the quan- 
tum Rabi model to characterize two coupling regions that we 
may call lower coupling region and higher coupling region, 
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respectively. The first region comprises the JC regime, where 
the RWA holds, together with the perturbative USC regime, 
where small deviations from the JC occur (g/u < 0.1). The 
second identified region comprises a jump towards the higher- 
coupling regime, g/u > 0.4, and forms a precursor of the 
DSC regime (g/u > 1). The intermediate region, that is for 
0.1 /u < g/u < 0.4, determines a kind of dark zone where 
no intuitive physics has been identified up to now. 

We begin by studying a suitable dynamical observable, the 
time average photon number 

1 /"* 

n = lim - / dt'n (t'), (2) 

which measures the break of the U (l)-symmetry in the JCM. 
As we will show, this quantity exhibits significant features to 
characterize the regions mentioned before. 

Since Hr commutes with the parity operator P, it is useful 
to investigate the dynamics within a fixed parity chain. Let \a) 
with a = ±1 correspond to the states \e) and \g) of the two- 
level system, respectively. Then the subspace H± with parity 
±1 is spanned by states {\<j> 8 ) ® | ± cr), \<j> a ) (g) | =p a)} where 
I ^(a), s) denotes the (anti-)symmetric part of |0) which is an 
element of the Hilbert space %b of the radiation mode. Hb is 
spanned by Fock states {|n)}, which are (anti-) symmetric if n 
is (odd) even. 

There exists a transformation F± that maps the element | 0) 
of 7-Lb onto a parity eigenstate |0, ±) that belongs to a parity 
chain 
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FIG. 1 . Panel (a) shows the time averaged photon number expecta- 
tion value no - the JC result of Eq. (6| and the exact result. Panel (b) 
shows the difference n§ xact — no°. The black line shows where the 
JCM matches exactly the analytical solution. 



|0, ±) = S ® | ± a) + a ® I T (3) 



The dynamical quantities in each chain depend only on the 
initial distribution of photons, which fixes |0(£),=b) at time 
t = 0. The state at later times follows as solution of the 
Schrodinger equation 



idM(t)) = H±m)), 



(4) 



where H± acts on functions in Hb B21 . The natural observ- 
able within the invariant subspaces is thus the photon number. 
We consider first the time-dependent expectation value, 

n </>W = (0(0 j -h|a^a|0(t), +) for some initial |0(O)) and 



positive parity. For A = 

|#0)> = \m) = (at)"7v 
at time t reads (g = g/o;) 



2 /2 = and an initial Fock state 
ra!|0), the average photon number 



n r 



l (t) = m + 2# 2 (1 - cos(cjt)), 



(5) 



which entails that the difference n m (t) — m is independent of 
the initial photon number, always greater than zero and oscil- 
lates with the mode frequency u. Because Fock states have 
definite reflection symmetry, the initial state in the full Hilbert 
space will be the product \m) <S> |(— l) m ) G H+. 

Starting now from the photon vacuum |0) we consider val- 
ues A ^ 0. Then, the correlation function no(t) depends on 
the initial state and shows a complicated oscillatory behavior. 
Still the time evolution can be globally characterized by the 
time average quantity defined in Eq. ([2]). Its dependence on 
A = A/u for fixed g shows features that allow to discern 
a regime which corresponds to the validity of the RWA (the 
JC-regime) from a region where the non-conservation of the 



excitation number in the QRM becomes relevant. As shown 
in Fig. [TJ a )> the Jaynes-Cummings model exhibits a sharp 
peak of no at A = 1/2 with vanishing width as g — > 0, 
and no < 1/2. The maximum value 1/2 is reached at reso- 
nance and corresponds to Rabi oscillations within the first JC- 
doublet, { |0e) , \ lg)}. One photon is periodically exchanged 
between qubit and radiation field which contains 1/2 photon 
on average. For general values of g and A, the value of no in 
the JCM reads 



2g 2 d 2 



(5 2 +d 2 ) 2 



(6) 



where d = A - | 



(A-f) 2 - 



On the other hand, the exact result exhibits an unbounded no 
due to the broken U (l)-symmetry, together with a pronounced 
broadening of the resonance. Figure [TJb) shows the differ- 
ence of the time-averaged photon number no calculated from 
the exact analytical solution and from the JCM. Notice that 
for A = 1/2, there is good agreement between both quan- 
tities until a value g ~ 0.1, establishing what we called be- 
fore the lower-coupling region. This means that with regard 
to photon production, the RWA result is valid for couplings 
associated with the USC, at least exactly at A = 1/2. For 
A 7^ 1/2 the deviations set in much earlier. Interestingly, 
for A > 1/2, less photons are generated than expected from 
the RWA calculation. This corresponds to a shift of the reso- 
nant A to values below u/2. Finally, for couplings g > 0.4, 
what we called before the higher-coupling region, the reso- 
nance is completely gone and the photon production exceeds 
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FIG. 2. In (a) we show the revival probability Po(t) of state 
1 0(0)) = |0) |e) and in (b) the revival probability P 4 (t) of state 
|0(O)) = |4) |e) for A = 0.25 and three values of g. The result 
depicted is obtained by solving Eq. (|4| for positive parity. 



the RWA prediction for all A, growing rapidly with the cou- 
pling. This marks the crossover to the deep strong coupling 
regime. Within the DSC, g > 1, the qualitative behavior of 
the system changes again, accompanied by the stabilization 



of Schrodinger cat-like states (see section IIIC). 



III. DEEP STRONG COUPLING FEATURES 

We study the time evolution of the autocorrelation function 
(0(O),+|0(£),+) where we start with an initial Fock state 
with 1^(0)) = \n). From this the revival probability is ob- 
tained as P n (i) = \(n, +|0(t), -f )| 2 . In Fig.[2j we observe the 
characteristic collapse and revivals of P n (t) in the deep strong 
coupling limit as in l26l . This feature is not present in the 
USC but appears only for sufficiently strong g. Starting in the 
vacuum as in Fig. |5Ja), the revivals are explained by the fact 
that the initial state |0) is a coherent state expressed in terms of 
the ground state of the quantum Rabi model |0) = |^o) 
for A = 0. For large g/A the physical behavior of the 
A = case is already present. That the physics of this limit 
is present already in a considerable distance from it is dis- 
cussed in Sec. |III A| Collapse and revivals are also observed 
for higher initial Fock states as seen in Fig.[2jb). Then the re- 
vivals become sharper but have additional contributions from 
higher harmonics. 

In Fig. [3] we show the time evolution of the distribution ob- 
tained by projection on the complete Fock state basis, i.e. the 
photon number distribution. Again the collapse and revival 
oscillations are visible in form of the red nodes of the interfer- 
ence patterns particularly pronounced in Fig.[3jd). The second 
feature that is observed is the "bouncing of photon number 
wave packets" l26l . Again the phenomenon can be under- 
stood in terms of the limit A = where the eigen basis of the 
quantum Rabi model becomes the basis of a shifted oscillator. 
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FIG. 3. The figure shows the time-dependent photon number distri- 
bution \{n\(j)(t))\ 2 for g = 0.7 (a,b) and g = 2.0 (c,d) and different 
initial Fock states |0(O)) = |0) <g> \e) (a,c) and |0(O)) = |16) <g> \e) 
(b,d). A = 0.25. The boundaries for the "bouncing of photon num- 
ber wave packets" can be obtained in a simple way using Fig.[4|b). 



This is discussed in detail in the next section, Sec. |III A] 



A. Physics of the adiabatic approximation 

The natural starting point for an approximation of the full 
quantum Rabi model at large g is the well known adiabatic 
approximation [33] . The Hamiltonian to zeroth order for fixed 
parity reads 



Hq = ua)a + g(a) 



(7) 



The characteristic Rabi physics in the DSC can already be ex- 
tracted by use of the eigenbasis { \n; g)} of this simple Hamil- 
tonian. Any more sophisticated approximation to the eigenba- 
sis improves the quality only marginally and does not reveal 
new physical behavior. For convenience we restrict ourselves 
to the positive parity subspace. The spectrum in is given 
by the set of zeros x n of the function [29], 



G+{x)= J2 R ™W 

ra=0 



rnuj 



(8) 
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as E n = x n — gg. The functional form of G+(x) reads 

g + (x) = g° + (x) + y 



(9) 



with G\(x) ~ (1 — A/x)e~^ and the h^(A) vanish for 
A = 0. G+(x) is of slow variation on the scale given by u, 
therefore the location of the zeros of G + (x) is determined by 

the poles at x = nuj^n = 0,1,2 This leads to an almost 

equidistant distribution of the x n for small A/g: The n-th root 
x n lies on the left or on the right of noo according to the sign 
of fo+(A). Because the h+(A) grow with A, we conclude 
that x n+ i — x n w uj for small A/g. The dynamics of the 
model for fixed parity is encoded in this smooth distribution 
of eigenvalues together with the fidelity of the adiabatic basis. 

There is a natural association of the shifted oscillator ba- 
sis {|n; g)} to the true eigenbasis that can be inferred 
from the representation [29], 



(10) 



oo A\fm\ 
\ipn) = ^2 K m(x n )- —\m;g). 



rn=0 



Xn — moo 



As the quantum Rabi model approaches the adiabatic limit, 
the spectrum converges to the spectrum of the shifted oscilla- 
tor, i.e. 



noo for g/A 



oo 



(11) 



This limit is encoded in the representation ( fT0| ) of the Rabi 
eigenstates by the divergence of 



(x n — noo) 1 — ^ oo for g/A 



oo. 



(12) 



This factor selects the nth shifted oscillator in the sum in 
Eq. ( fT0| ) for sufficiently large g/A. We conclude that the Rabi 
eigenbasis in a fixed parity subspace is close to diagonal in the 
basis of shifted harmonic oscillators. 

This behavior is illustrated in Fig. |4ja) where we plot the 
contributions from the projections of the eigenstates on the 
shifted oscillators \(m; g\i(j+)\ 2 on a logarithmic scale. The 
resulting matrix shows a rapid decrease of the off-diagonal 
elements. This reduction is already present for small values 
of n and m and becomes more pronounced for higher values. 

In Fig. |4jb), by contrast, we plot the representation of 
the eigenstates in terms of Fock states. The figure dis- 
plays a clear parabolic shape. This shape can be com- 
pletely understood within the adiabatic limit. Considering 
the representation of the shifted oscillator states in terms of 
Fock states, a perfectly symmetric parabola is seen. This is 
due to the fact that the Fock basis is transformed into the 
shifted oscillator basis by the operator D(—g) = e~^ a 
and the group property of D, D(x)~ 1 = D{—x). Be- 
cause the matrix element D mn {— x) reads (m;0\x;n) = 
e- x2 ' 2 x m - n ^m\/n\L 7 ^ rn {x 2 ), where L^ m denotes a La- 
guerre Polynomial, it follows \D nm (x)\ = \D nm (—x)\, fur- 
thermore, D nm (-x) = D~^(x) = D*(x)m n and there- 
fore \D nm (x)\ = \D rnn (x)\. The finite value of A = 0.25 
then leads to the slight asymmetry visible in the parabola of 
Fig. |4p}). Surprisingly, also Fig. |4ja) shows an almost sym- 
metric shape which cannot be explained by the preceding ar- 
gument. 
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FIG. 4. Density plot of the projection of eigenstates onto shifted 
oscillator states | (m; g\ip n ) | 2 (a) and Fock states | (m; 0\ip n ) | 2 (b) 
for g = 0.7 and A = 0.25. This is for positive parity, i.e. ip n is an 
eigenstate of H+. The arrows in panel (b) show how to determine 
the boundaries for the "bouncing photon number wave packets" that 
are visible in Fig. [3] 



The effect of bouncing photon number wave packets men- 
tioned earlier and shown in Fig. [3] can be described solely by 
the presence of the parabola in Fig.|4jb). To illustrate the de- 
termination of the boundaries for the bouncing phenomenon 
in Fig.[3p3), we have added the arrows in Fig.|4jb). The up- 
wards pointing arrow indicates the initial state |16). The hori- 
zontal arrows then lead to the minimum and maximum photon 
numbers contained in the eigenstates. 

In Fig. |5j we present results for the time evolution of the 
quantum Rabi model calculated using the adiabatic approxi- 
mation. Fig.^b) shows that for g/A = 8 the adiabatic ap- 
proximation gives a detailed picture of the full time depen- 
dence. Fig. [5ja) by contrast shows that for moderate values 
of the coupling, here g/A ~ 2.8, the full adiabatic approx- 
imation, i.e. approximated spectrum and approximated ba- 
sis, completely misses a physical description of the dynam- 
ics. This insufficiency of the full adiabatic approximation is 
mainly due to deviations in the spectrum which yield wrong 
phases. By contrast, the approximated basis combined with 
the exact spectrum leads to a qualitatively correct description 
of the dynamics. 
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FIG. 5. The time evolution of the photon number expectation 
value rio(£) for g = 0.7 (a) and g = 2.0 (b). no(t) is eval- 
uated as sum over the Photon distribution of Fig. [3ja) and (c): 
M*) = J2m m \(Mt)\m)\ 2 where |0 o (t)) - e - <H +*|0> = 



y 



i |^n)(V ; n|0). (i) solid lines depict numerically exact re- 



sults calculated with this expression, (ii) dashed lines are calculated 
with \ip n ) replaced by \n,g) and (iii) dotted lines are calculated with 
the additional replacement of the spectrum by its first order corrected 
perturbation theory value used in Eq. (14) . For g = 0.7, the adiabatic 
approximation allows a good description of the time evolution only if 
the exact spectrum is employed as seen in panel (a). For g = 2.0, the 
adiabatic approximation is highly precise also if the approximated 
spectrum is used. 



To further illustrate the point that the adiabatic basis con- 
tains already all of the physics even far away from the limit 
g/A — )> oo, while the adiabatic spectrum leads to unphysical 
results, we define a distance of the state |n, g) of the approxi- 
mate basis from the true eigenstate \i/j n ) by 



D n (g,A) = l-\(n;g\^ n )\ 2 



(13) 



A distance of the first order corrected eigen energy from the 
true eigen energy is given by the difference 

D*(g, A) = \un - gg + A(-l) n (n, g\n, -g) - E n \. (14) 

These functions are shown in Fig. [6] for n = 0. While their 
magnitude cannot be compared directly, their dependence on 
g and A allows to conclude the following. The slope of the 
contour lines is much steeper in case of the basis (Fig. [6ja)) 
compared to the spectrum (Fig. |6jb)). This means that the 
quality of the approximate basis increases much faster with 
g/A than that of the first order spectrum. This is surprising, 
as the basis is only a zeroth order approximation while the 
spectrum is approximated to first order. 

It is possible to improve the adiabatic basis in a systematic 
fashion (34j [33 or use the RWA on top of the adiabatic ap- 
proximation (36). There are other approaches which reduce 
the approximative diagonalization of Hr mH± to diagonal- 
ization of finite matrices |37] [38). All these techniques are 
equivalent to a break-up of H± into a set of invariant sub- 
spaces with finite dimension. This means that a continuous 
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FIG. 6. The distance A)(<7, A) of the shifted oscillator ground state 
from the ground state of the quantum Rabi model in the positive par- 
ity subspace is shown in panel (a). The corresponding difference of 
the eigen energies Do(g, A) in panel (b). The contour lines in the 
plots are drawn at values of 0.01 and 0.05. Comparing the slope of 
the contour lines between panel (a) and (b) shows that the quality of 
the basis increases faster than the quality of the spectrum. 



symmetry is superimposed on the quantum Rabi Hamilto- 
nian, leading to an effective model of Jaynes-Cummings type 
|29) . In the quantum Rabi model, such an additional (hidden) 
symmetry is not even present in an approximate sense, be- 
cause otherwise Fig.[4ja) would display a block-diagonal pat- 
tern. Instead, the off-diagonal components of the true eigen- 
states in terms of shifted oscillator states are distributed rather 
smoothly on both sides of the central diagonal being at the 
same time smaller by at least two orders of magnitude. It 
can be concluded that no self-consistent reduction to finite- 
dimensional invariant subspaces (36J [33 is possible which 
improves the adiabatic basis without creating a strong (and 
unphysical) symmetry. Most of the physics in the DSC regime 
(in fact, already for g > 0.7) is captured by the simple adia- 
batic basis. 

Besides these physical considerations, we point out the 
technical fact that the spectrum is much easier to compute than 
the basis. Combining the exact spectrum with the simple ba- 
sis allows us to predict correctly dynamical quantities of the 
QRM in the blue sketched parameter region of Fig. [6ja), as 
demonstrated in Fig. [5] 
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FIG. 7. Wigner functions of eigenstates in the parity chain 
and considering the DSC regime, (a) Grouns state, (b) first excited, 
(c) second excited, and (d) third excited states. We can see that they 
resemble the Wigner functions of Fock states |0), |1), |2), |3), re- 
spectively, but displaced from the center as expected from the shifted 
oscillator basis. 



B. Eigenstate Wigner functions 
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FIG. 8. Revival probability P a (t) of coherent states e aa ^ |0) for 
different values of a and A = 0.25. If a = g, the time evolution of 
the coherent state remains almost constant, P a =g (t) = Pc + (t) ~ 1 
for all times. This is to be compared with the time evolution depicted 
in Fig.JTJb) of the initial Fock state |4) that yields the same photon 

number (h) = 4 as \C + ) = e 2at |0). 



The above discussion in the DSC regime led us to find the 
shifted basis as a good approximation to the real one. The fea- 
ture can be further confirmed by analyzing the Wigner func- 
tions associated to each eigenstate in both parity chains, T~L±. 
Fig. [7] shows the Wigner functions corresponding to the first 
eigenstates in the parity chain We see that these func- 
tions resemble the Fock states |0), |1), |2), and |3), though 
shifted from the center, as expected from the shifted oscillator 
basis. Notice the appearance of some additional interference 
pattern coming from other Fock states which contribute to the 
whole solution, see Fig.[4ja). It is noteworthy to mention that 
exactly the same behavior can be observed in the other par- 
ity chain %+. Both chains behave very similar regarding the 
form of the eigenvectors, the ground state in %_ looks almost 
the same as the first exited state in %+. 



C. Stationary Schrodinger cat-like states 



ton number (n) in e ± ^ a |0) is g 2 , an experimentally realiz- 
able value of (n) = 4 corresponds to a DSC value g — 2. 
These states behave almost like eigenstates regarding expecta- 
tion values of observables, e.g. the revival probability Pc± (t). 
As seen in Fig. [8] Pc± (t) ~ 1 remains almost constant for 
a = g and exhibits only small high-frequency variations, but 
not the collapse-and-revival behavior with frequency uj as the 
Fock states in Fig.[T] 

Fig. [8jb) shows that, for a = g, the oscillations with u 
are completely gone and all fluctuations occur on time scales 
much shorter than the oscillator period 2tt/uj. Averaged over 
such a short time, Pc± (t) remains constant and behaves effec- 
tively a conserved quantity although the corresponding state is 
not an eigenstate. The phase variable of the protected coher- 
ent states depends directly on the coupling g and confines the 
effect to the DSC regime if the coherent states are required to 
contain more than two photons on average. 



Fig. |6ja) shows that the ground state of H+ at strong cou- 
pling is well approximated by |0; g) = V 2 -^ 1 " |o). If the 
system is prepared in this state, it will only weakly depend 
on time. Using the transformation F± from Eq. ([5J, the state 
|0; g) is mapped to the following states in the original basis 
with fixed parity ±1: 

\C + ) = e-^ 2 / 2 (cosh(^at)|0) \e) - sinh(£at)|0) \g)) 
|C_) = e-^ 2 / 2 (cosh(^at)|0) \g) - sinh(£at)|0) \e)) 

(15) 

These Schrodinger cat-like states contain qubit-field entangle- 
ment (39), where the terms sin h(gd*)\0) and cosh(^at)|0) 
are (anti-)symmetric superpositions of the semiclassical states 
e-^ at |0) and e^ at |0) B51BH . Indeed, as the average pho- 



IV. CONCLUSIONS 

We have studied the dynamical properties of the quan- 
tum Rabi model. By using the recently developed analyt- 
ical solutions 1 29], we obtained qualitative features of the 
spectrum and the eigenstates in two parts of the QRM: a 
lower and a higher coupling region, respectively. It turns out 
that it is mandatory to pay attention to the sole and impor- 
tant Z2 -symmetry of the QRM, which leads to the separation 
of the Hilbert space into two invariant subspaces, the parity 
chains l26l . Within each chain, the dynamics appears to be 
rather simple. Especially in the deep strong coupling regime, 
g > 1, we find a quite regular periodic behavior of the photon 
number distribution which can be traced back to the almost 
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equidistant separation of energy levels. This is no longer true 
if the initial state has no definite parity. Under such circum- 
stances, e.g. if the initial state is a product of a coherent state 
and a qubit eigenstate, the time development may be more 
complicated, as both chains interfere l32l . 

We have identified a simple observable, the average photon 
generation from the vacuum at positive parity, to discern the 
convenient lower-coupling and higher-coupling regions. This 
quantity shows a remarkable sharp peak just for A = lj/2, 
which demonstrates in a direct way the resonant enhancement 
of coupling between qubit and cavity mode for small g, which 
was the original motivation to introduce the RWA in the first 
place El. This determines the lower-coupling region. For 
stronger couplings, the resonance gets broader until the peak 
vanishes around g/uj ~ 0.4, establishing the higher-coupling 
region. Here, the average generated photon number becomes 
larger than the JC limit of 1/2 due to the counter-rotating terms 
in ([]]), which break the conservation law of the JCM. 

As early as for g/uj > 0.7, the characteristic features of the 
deep strong coupling regime begin to manifest. The dynamics 
for fixed parity are dominated by the adiabatic basis and the 
almost equidistant spectrum. Interestingly, the nontrivial ef- 
fects, which separate the QRM in this regime from the simple 



adiabatic limit, can be incorporated by using the exact spec- 
trum together with the adiabatic basis. In fact, the exact eigen- 
states are very close to their adiabatic approximants whereas 
deviations in the eigenvalues lead to phase differences which 
become apparent after longer times (Fig. [5]). The fidelity of the 
adiabatic basis is even more visible in the Wigner representa- 
tion, which demonstrates the similarity of parity eigenstates 
with Fock states in the deep strong coupling regime (Fig. [7]). 

Finally, the DSC allows for a special class of states which 
are not exact eigenstates but lead to expectation values fluctu- 
ating with small amplitudes on very short time scales. In this 
sense, they form stationary Schrodinger cat-like states, unaf- 
fected by the system interaction. 
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